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We calculate, in the framework of QCD sum rules and to next-to-leading order in 
perturbation theory, the universal function £3(1; • v') which appears at order 1/mg in 
the heavy quark expansion of meson weak decay form factors. We find that radiative 
corrections of order a s are very important. Over the kinematic range accessible in 
semileptonic decays, ^(v-v 1 ) is proportional to the leading-order Isgur-Wise function 
£(v ■ v') to very good accuracy. Taking into account all sources of uncertainty, we 
estimate £3/*$; = (0.6 ±0.2). This reduces the theoretical uncertainty in the extraction 
of I Vcb\ from B — > D Iv transitions. A measurement of the form factor ratio A2/A1 
in B — > D*£ v decays can be used to test our prediction. 
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I. INTRODUCTION 



In the limit where the charm and bottom quarks are considered infinitely heavy, their 
strong interactions with light quarks and gluons acquire additional symmetries [0-^]. In 
particular, the weak decay form factors describing the semileptonic transitions B — > D Iv 
and B — > D*£v become related to a single universal function. The normalization of this 
so-called Isgur-Wise form factor is known in the zero recoil limit, where the initial and 
final meson have the same velocity. This allows a model-independent determination of 
the Cabibbo-Kobayashi-Maskawa matrix element | V c \\ up to corrections arising from 

the fact that m c and m& are, after all, not infinitely heavy. In this work, we investigate 
a particular type of such corrections, which is important to the extraction of | V c b\ from 
B — > D iv transitions. 

The heavy quark effective theory (HQET) provides a convenient framework to analyze the 
weak decays of hadrons containing a heavy quark [P-[T2]| . It provides a systematic expansion 
of hadronic matrix elements in powers of I/itlq. The coefficients in this expansion are 
universal functions of the velocities v and v' of the initial and final hadrons (v 2 = v' 2 = 1, 
v ■ v' > 1), but they do not depend on the flavor or spin of the heavy quarks. These form 
factors originate from long-distance hadronic dynamics, so they can only be investigated 
using some nonperturbative approach to QCD. One such method is provided by QCD sum 
rules ||13|| , which have recently been widely used to calculate hadronic matrix elements in 
HQET |l4|-p2||. At leading order in the heavy quark expansion, a single Isgur-Wise function 
^{y-v',n) is required to parameterize the current-induced transitions M{y) — > M'(V), where 
M and M' are pseudoscalar or vector mesons containing a single heavy quark This is 
conveniently expressed by the compact trace formula 

(M'(v')\h'Th\M(v)} = -^(v-v'^)Tt{M'{v')TM(v)}, (1) 

where T is an arbitrary Dirac matrix, and h and h! denote the velocity-dependent effective 
fields in HQET which represent heavy quarks Q and Q' moving at the hadron's velocities v 
and v' . The heavy mesons are represented by covariant tensor wave functions 

, N . (1 + 0) f — 75 ; pseudoscalar meson, 

MM = Jtum- — { j , (2) 

v 2 I f ; vector meson, 

which have the correct transformation properties under Lorentz boosts and heavy quark spin 
rotations. Here myi denotes the physical meson mass, and e is the polarization vector of 
the vector meson. Current conservation implies that the Isgur-Wise function is normalized 
at zero recoil: £(l,/x) = 1. Except at this point, the universal form factor depends on 
a subtraction scale /i, since the velocity-changing currents in the effective theory have to 
be renormalized. The /^-dependence of the Isgur-Wise function cancels against that of the 
Wilson coefficients which appear in the matching of currents of the full theory onto currents 



of the effective theory 



At order l/m,Q, matrix elements receive contributions from higher dimension operators 
in the effective Lagrangian and in the effective currents |Tl| . The former give rise to three 
new universal functions, usually denoted by Xi( v ' v ' \ A 4 ) f° r i — 1)2,3. In the framework 



of QCD sum rules, we have investigated these form factors in Refs. [19 , 21| - Here we shall 



focus on the second type of corrections, which come from operators of dimension four in 
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the short-distance expansion of the weak currents in HQET. In the case of the vector or 
axial vector currents, there are fourteen independent operators of this type. Their Wilson 
coefficients have recently been calculated to next-to-leading order in renormalizat ion-group 
improved perturbation theory p5 |. On dimensional grounds, any dimension four current 



operator must contain a covariant derivative acting on one of the heavy quark fields. Thus, 
these operators have the generic form h' V iD a h or (iD a h') V h, where T is again an arbitrary 
Dirac matrix. In analogy to (|T|), the corresponding matrix elements can be parameterized 
by a tensor form factor £ a defined by pT] 



(M'(v')\tiTiD a h\M(v)) = -ATi{^ Q (v,v',fi)M'(v')TM(v)}, (3) 

where A = mu — rriQ = tum' — vtlqi denotes the mass difference between a heavy meson 
and the heavy quark that it contains. This parameter sets the canonical scale for power 
corrections in HQET. Matrix elements with a derivative acting on h! are related to (|3|) by 
Dirac conjugation. The most general decomposition of the form factor can be written 

£, a (v, v', n) =i+{y ■ v', fi) + v') a + £-(v- v', p) (v - v') a - £ 3 (v • v', fi) >y a . (4) 

Due to T-invariance of the strong interactions the coefficient functions £j(t> • v',fi) are real. 
Furthermore, the equation of motion of HQET can be used to derive the constraints JTT[ 



I v - v ' - I 1 
£ + (v ■ v, fi) = - v . v , + 1 • v', fi) - v . v , + l Uv ■ v', fi) , (5) 

These relations show that only one of the coefficient functions, say £a(v-v', fi), is independent. 
They also suggest a close relation between the subleading functions £j(t> ■ v', //) and the 
leading-order Isgur-Wise function £(v ■ v',fi). The origin of this relation is the invariance 
of the effective theory under reparameterizations of the heavy quark momentum pp] , pB| . In 
fact, one can show that the /i-dependence of the functions £j(t> • v',fi) is the same as the 
//-dependence of the Isgur-Wise function. This leads us to introduce a new, renormalizat ion- 
group invariant function i](v ■ v') by 

V{v-v) = — r — . 6 

We expect that r)(v ■ v') is a slowly varying function of order unity. By means of (|5|) and 
(§), meson matrix elements of the dimension four operators in the short-distance expansion 
of the currents can be parameterized in terms of the product A£(v ■ v',fi) and a single new 
function i](v ■ v'), which is independent of the renormalization scale. 

Knowledge of the function r\ (v ■ v') becomes important when one wants to extract | V c b\ 
from a measurement of the differential decay rate dT/d(v ■ v') for B — > D Iv transitions near 
zero recoil. Because of the known normalization of the Isgur-Wise function at v ■ v' — 1, 
hadronic uncertainties affect such a measurement only at the level of power corrections of 
order l/m c and l/m^. Note that whereas the B — > D*lv decay rate is protected against 



first-order power corrections by Luke's theorem [11], the rate for B —> D iv is not, due to its 
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helicity suppression at zero recoil p7| . However, one can show that the symmetry-breaking 
corrections are parametrically suppressed by the "Voloshin-Shifman factor" || 



S 



( m B - m D \ 2 
\m B + m D J 



0.23 



(7) 



For this reason one may hope that the theoretical uncertainty in extracting | V c b\ from these 
transitions is not much worse than in the case of B — > D*iv decays An extrapolation 
of the spectrum to zero recoil gives 



lim 



dT(B -> D£P) 



Gl I V cb \ 2 



v-v — >1 



where 



K = 5 1 + 



d(t> ■ 


v') 


f A 


+ A 


\2m c 


2m,bJ 



48tt 3 



jn B + m D ) 2 m 3 D 



l + S-K 



(8) 



(l + 5i)- 2 (1 + $3)77(1) 



(9) 



Here 5i are perturbative QCD corrections arising from finite renormalizations of the currents 
in the intermediate region mj, > > m c . Numerically, one finds 5i ~ 6% and 5 2 w 9% 
24] , |2o| . The only unknown quantity (except the heavy quark masses) in this equation is the 
value of the function rj(v ■ v') at zero recoil. 

It has been pointed out in Ref. [IT]] that certain ratios of the hadronic form factors 
describing B — > D*iv transitions are very sensitive to symmetry-breaking corrections to 
the heavy quark limit. Their measurement can be used to test model calculations of the 
subleading universal form factors. Consider the quantity 



1 - 



AM 2 ) 



(10) 



where A\(q 2 ) and A 2 {q 2 ) are axial vector form factors in the notation of Ref. [p8 |. Introducing 
the variable v-v' instead of q 2 , and performing a 1/ itlq expansion of the hadronic form factors, 



one obtains [19 



1 - 



A 



I 3 \ 

, , , — + — )r](v -v') + ... . 
v • v + 1 \m c nib J 



(11) 



The ellipsis represents a small "hyperfine correction" proportional to the function X2{v-v', fi), 
and higher-order power corrections of order I/ttlq. The perturbative corrections to R 2 turn 
out to be completely negligible. Unless the function r](v ■ v') were highly suppressed, the 
dominant symmetry-breaking correction to R 2 is the one shown in ([TT|) . A measurement of 
this ratio would therefore provide valuable information about the form factor r](v ■ v'). 

the function £ 3 (v-v', /i) has been analyzed using the QCD sum rule approach 



In Ref. 19 



and adopting the standard approximations, in which one neglects radiative corrections. Such 
an approach leads to the parameter-free prediction that £3(1, /i) = r/(l) = 1/3. In the context 
of QCD sum rules, corrections to this simple result can only come from radiative corrections 
or higher dimension condensates. It was recently found that such corrections can be quite 
significant, however. For the subleading Isgur-Wise functions x^(v ■ v',fi) and X3( v ' ^'jAOj 
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the terms of order a s even give the dominant contributions to the sum rules In the 
case of £ 3 (v ■ v',fi) this may be expected to be even more so, as due to the peculiar trace 
structure associated with this form factor in (^) there is no quark condensate contribution at 
tree level. Therefore, it is important to refine the sum rule analysis of Ref. Jl9| by including 
radiative corrections. This requires, in particular, the calculation of the two-loop corrections 
to the triangle quark-loop diagram. The techniques to handle the corresponding two-loop 



integrals have been developed in Ref. 



II. DERIVATION OF THE SUM RULE 

The QCD sum rule analysis of £ 3 (t> ■ v', //) proceeds in complete analogy to that of the 
Isgur-Wise function. For a detailed description of the procedure the reader is referred to 
Refs. [jl9p2"l . Here we shall only very briefly sketch the main steps. One considers, in the 
effective theory, the three-current correlator 

J dx dx' e ^'-'- fc -) ( | T { [qT M ,h%„ [h' T iD a h} , [h T M qQ | ) 

= Tr{ Ea(v, v', k, k') Tw^-T ^ T M } , (12) 

where v and v' are the velocities of the heavy quarks, and k and k! are the external off-shell 
momenta injected into the three-point function. Depending on the choice Tm = —75 or Tm = 
Ifi ~ v fi, the heavy-light currents interpolate pseudoscalar or vector mesons, respectively. 
The Dirac structure of the correlator, as shown in the second line, is a consequence of heavy 
quark symmetry, as reflected in the Feynman rules of HQET. The quantity H Q obeys a 
decomposition analogous to with coefficient functions S± and H 3 that are analytic in 
the "off-shell energies" 00 = 2v ■ k and uj' = 2v' ■ k', with discontinuities for positive values of 
these variables. These functions also depend on the velocity transfer y = v ■ v'. From now 
on we shall focus on the coefficient E 3 (u,uj',y), which is used to construct the sum rule for 
the subleading form factor £3(2/, /i). 

The idea of QCD sum rules is to relate a theoretical approximation to the operator 
product expansion of S3 to a hadronic representation of the correlator in terms of physical 
intermediate states. The lowest-lying states are the ground-state mesons M and M' asso- 
ciated with the heavy-light currents. They lead to a double pole located at uj = uj' — 2A. 
The residue is proportional to the function £3(2/, /x). One finds ||T9| 



^poic/ / x A £3(1/,//) F 2 (/i) 

(u — 2 A + ie)(uj' — 2A + it) 

where F corresponds to the scaled meson decay constant (F ~ /m\/™m)' Both F and 
£3 are defined in terms of matrix elements in the effective theory and are scale-dependent 
quantities. In the deep Euclidean region, the correlator can be calculated perturbatively be- 
cause of asymptotic freedom. The idea of Shifman, Vainshtein, and Zakharov was that, at 
the transition from the perturbative to the nonperturbative regime, confinement effects can 
be accounted for by including the leading power corrections in an operator product expan- 
sion. They are proportional to vacuum expectation values of local quark-gluon operators, 
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the so-called condensates fL3| . Following the standard procedure, we write the theoretical 
expression for H 3 as a double dispersion integral and perform a Borel transformation in uj 
and uj'. This yields an exponential damping factor in the dispersion integral and eliminates 
possible subtraction polynomials. Because of the flavor symmetry of HQET, it is natural 
to set the associated Borel parameters equal: r = r' = 2T. Following Ref. ||20|| , we then 
introduce new variables uj + = \ (u + uj') and uj- = uj — uj', perform the integral over u)-, and 
employ quark-hadron duality to equate the remaining integral over uj + up to a threshold ujq 



to the Borel transform of the pole contribution in (0). This yields the Borel sum rule 



^0 

A UV, AO F 2 (p) e~ 2K ' T = J duj + e~^ T p(u + , y) = K(T, uj , y) 



(14) 



The effective spectral density p arises after integration of the double spectral density over 
uj-. For practical purposes it is useful to notice that the uj + - dependence of each term in 
p(u>+,y) is known on dimensional grounds [21,22]. It thus suffices to calculate the Borel 
transform of S3, corresponding to the limit ujq — > 00 in flTil) . The dependence on ujq can be 
reintroduced later. 

As pointed out above, the theoretical expression for the right-hand side of the sum rule 
consists of a perturbative part and condensate contributions: K = A' pert -|-A' cond . Let us first 
present the result for the nonperturbative contributions. The leading terms are proportional 
to the quark condensate (dimension d — 3), the gluon condensate (d — 4), and the mixed 
quark-gluon condensate (d — 5). For a consistent calculation at order a s , we calculate the 
Wilson coefficients of the quark and gluon condensates to one-loop order, and the coefficient 
of the mixed condensate at tree level. The truncation of the series of power corrections at 
the mixed condensate seems safe. The contributions from four-quark operators (d = 6) are 
suppressed relative to the quark condensate by a factor \ (qq)\/T 3 ~ 1 — 5%. The calculation 
is most conveniently performed using the coordinate gauge x ■ A(x) = with the origin 
chosen at the position of the velocity-changing heavy quark current. We find 



K cond (T, 00, y) 



where 



2a s (qq)T 
3tt 



r(y) 



2 + r{y) 



+ 



a s GG)/y-l\ (qg s VapG a(3 q) 



96vr 



y + 1 



12T 



in (y + yy 2 - i 



(15) 



(16) 



Let us now turn to the perturbative contributions to the sum rule. At order a s , one has 
to evaluate the bare quark loop as well as the seven two-loop diagrams depicted in Fig. [1|. 
A new feature of the present sum rule, as compared to the sum rule for the Isgur-Wise 
function considered in Ref. p2| , is that there is a diagram (-D7) where a gluon originates 
from the covariant derivative contained in the current. We denote the Borel transformed 
contributions of the individual diagrams to the function K peit (T, 00, y) by L>i. Throughout 
the calculation we use Feynman gauge and dimensional regularization. The bare quark loop 
is readily calculated and gives |19| 



Do 



2N C T D 
(4ti) d / 2 Ky 



D/2 



(17) 
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where D is dimension of space-time. The evaluation of the two-loop corrections is more 
complicated. For a detailed and systematic discussion of the techniques used to calculate 



the two- loop integrals the reader is referred to Ref. |22|. The present calculation proceeds 
very similar to that in this reference, where the sum rule for the Isgur-Wise function was 
derived at two-loop order. It is convenient to introduce a constant 



.4 



AN c C F g 2 s (2T) 



\2D-6 



(4tt)^ [ 2 (y + 1)] 



D-2 



-1) 



(18) 



where Cf = (N% — 1)/2N C . Once this quantity is factored out, we find that the contributions 
of the first four two-loop diagrams in Fig. [T] are the same as the corresponding contributions 
to the sum rule for the Isgur-Wise function: 



EA = 4 



- yr(y) 



+ 2 



yr(y) In [2(1 + y)] + 2yh(y) - 4 + 0(e) 



(19) 



We have expanded the result in e = (D — 4)/2 and introduced the function 



1 



L 2 (l - yi) - L 2 (l - y_)j + - y/y 2 - 1 r\y) 



(20) 



where y_ = y — \fy 2 — 1, and L 2 (x) is the dilogarithm. The calculation of the remain- 
ing three diagrams is more cumbersome. It requires rather elaborate techniques such as 



Kotikov's method of differential equations f29| . However, remarkable simplifications take 
place as we add up the various contributions. In particular all dilogarithms, which appear 
in intermediate steps of the calculation, cancel out. The final result is rather simple: 



2=5 



.4 



2tt 2 



(y 2 -l)r 2 (y)-(y + l) 2 + r(y) + 0(e) 



(21) 



Except for the last term, this is again the same result as for the Isgur-Wise function. 

The ultraviolet divergent terms in the sum of the seven two-loop diagrams match with 
the anomalous dimensions of the heavy-heavy and heavy-light currents contained in the 
three-current correlator in (|12"D. Thus the 1/e pole disappears upon renormalization of the 
currents. In the modified minimal subtraction scheme, the renormalization factors are |3,10 



'hi 



27re ' 



where 



11 , 4vr 

- = ~ + 1e - In —r 

e e n z 



(22) 



(23) 



and 



i(y) = 3 [yr(y) 



(24) 



Our exact two-loop result for the renormalized perturbative part of the correlator is 



7 



K peit (T, oo, y) = Z^Z U 2 i% + £ A 
3T 4 / 2 



i=i 



it 2 \y + 1 



1 + 



7T 



2-7(1/) ( ln^ + 7 £-y 



+ Cpert(3/) + Ml/) 



4tt 2 



17 

t 



(25) 



where 

Cpevt(y) 

Sc(y) 



i{y) 



y+l 



41n2-3 + ln 



- + - (y + 1) [2 + r(y) 



y%) - 1 



+ \n y -^ + \tf-l)r\y) 
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3/r(j/) 



(26) 



The function c per t(y) is the same that arises in the calculation of the sum rule for the Isgur- 
Wise function. 

The final expression for the QCD sum rule (|14"D is obtained when we reintroduce the 
continuum threshold co>o and write the result as a dispersion integral. This gives 



167T 2 \y + 1 



x < 1 H 

I 7T 



/ \1 , Z 1 47r 2 17 . . . , , 

2 - 7(y) In — + — + — + c peTt (y) + 6c{y) 
J co>_i_ y o 



2 + r(y) 



3tt 
12T 



dcj i e 



96tt \y+l 



(27) 



It is instructive to compare this to the sum rule for the product A£(y,/j,), which can be 



obtained from the two-loop calculation of Ref. |22 

2 ^ 



I67H \y + 1 



d^ e- w+/T uol 









x |l H 








71 





M i, // 4vr 2 17 . , 

2 - 7(y) In — + — + — + c pcvt (y) 
J uja_ y o 



2« 3 (gg) 
3vr 



yr(y) 



(a s GG) (y - 1 
96tt Vy + 1 



(2y + l) (qg s cr aP G^q) 



AT 



(28) 



Notice that the /i-dependence in (p7|) and ( p8|) is the same. Thus, our explicit calculation is 
in accordance with the fact that the function rj(y) in ([]) is renormalization-group invariant. 
We write 
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v(y) 



1 



and find that A(y) obeys the sum rule 

a,T 4 / 2 



+ A(y) 



A(y) Ae(z/)^ 2 e- 2A / T 



where 



12tt 3 Vy+ 1 
2a s (qq)T 



2 [ll + Qy+(3 + y)r(y)]5 3 (^ 



+ 



9vr 

(a s GG) f y-1 
y + 1 



7+(3-y)r(y) 



72tt 



(qg s <?ai3G al3 q) 



18T 



(y-i) 



(29) 



(30) 



5 n (x) 



r(n + i) 



dz z n e 



(31) 



Note that, since the right-hand side of (|30|) is of order a s , in this sum rule one is not sensitive 
to the running of the quantities £(?/) and F . Their /^-dependence would show up at order 
a 2 . For the analysis of the sum rule it is, therefore, consistent to use 



H(y)F 2 e 



-2A/T 



9T 4 



y + 1 



T 



AT 



(32) 



which is obtained from (|2q ) by neglecting terms of order a s . Taking the ratio of ( p0|) and 
reduces to a minimum the systematic uncertainties in the calculation of A(y). 



III. NUMERICAL ANALYSIS AND CONCLUSIONS 

In its final form, the sum rule for A(y) very much resembles the sum rules for the other 



subleading Isgur-Wise functions Xziv) and Xz{y)-> which we derived in Ref. [^]. Accordingly, 
the numerical analysis proceeds in a similar way. For the QCD parameters we take the 
standard values 

(qq) = -(0.23 GeV) 3 , 
(a s GG) = 0.04 GeV 4 , 
(qg.(Taf)G af> q) = m 2 (qq) , m 2 = 0.8 GeV 2 . (33) 

Our results turn out to be very stable against variations of these numbers within reasonable 
limits. Furthermore, we use a s /ir = 0.1 corresponding to a scale \x ~ 2A rs 1 GeV, which is 
appropriate for evaluating radiative corrections in the effective theory. Combining (|30|) and 
(|32|), we obtain A(y) and hence rj(y) as functions of Uq and T. These input parameters can 
be determined from the analysis of a QCD sum rule for the correlator of two heavy-light 
currents in the effective theory One finds good stability for u = 2.0 ± 0.3 GeV, 

and the consistency of the theoretical calculation requires that the Borel parameter be in 
the range 0.6 < T < 1.0 GeV. 
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In Fig. 0(a) we show the zero recoil value of the form factor i](y) as a function of the 
Borel parameter, for three different values of the continuum threshold. For y = 1 the sum 
rules simplify considerably. We find 



/u \ 8vr 2 (qq) /w \ 

Ay) ~~^ do \T ] 

3 ' 9tt 5 2tt 2 ml (qq) 



^ 1 ^ -W-T^W , (34) 



T y 9 T 5 



Neglecting the terms of order a s , we would recover the result 7/(1) = 1/3 derived in Ref. [19 
However, as seen in the figure these contributions are by no means negligible. They enhance 
the form factor by almost a factor 2. It supports the self-consistency of the sum rule approach 
that we find stability in essentially the same region of parameter space that leads to stability 



of the two-current sum rules considered in Refs. 17,181, and of other three-current sum rules 



analyzed in Refs. [21,22 



Over the kinematic range accessible in B — ► D^*>tv decays, we show in Fig. |](b) the 
range of predictions for rj(y) obtained for 1.7 < u < 2.3 GeV and 0.6 < T < 1.2 GeV. The 
numerical analysis confirms our guess that rj(y) should be a slowly varying function of order 
unity, which was the motivation for its introduction in the first place. In fact, the sum rule 
predicts that r)(y) ~ 0.6 essentially independent of y. The main uncertainty comes from 
the values of a s and a>o, which are not very accurately known. However, one should keep in 
mind that there are systematic uncertainties inherent in QCD sum rules which cannot be 
estimated by simply varying the input parameters. To be conservative, we quote our final 
result as 

r)(y) « 0.6 ± 0.2 ; 1.0<y<1.6. (35) 

This result has important implications for the extraction of | V c b\ from B —> D tv decays. 
According to fl9|), the I/tuq corrections to the decay rate are proportional to [1.06—2.18 7/(1)], 
and by a fortunate accident this combination is strongly suppressed for 7/(1) in the range 
(|35|). For the symmetry-breaking corrections to the decay rate in ([8]), we obtain (we use 
A = 0.5 GeV, m c = 1.45 GeV, and m b = 4.8 GeV) 

1 + S ■ K » 1.00 ± 0.03 , (36) 



i.e., at most a few percent. This is comparable to the expected size of 1/itlq corrections [|30| 



We conclude that the theoretical uncertainty in the determination of | V c b\ from this decay 
mode is not worse than in B —>■ D*iv transitions. Of course, the experimental measure- 
ment of B — » Div near zero recoil is more difficult. The reward of such a measurement, 
however, would be an independent determination of | V c b\ with surprisingly small theoretical 
uncertainties. 

Let us finally point out how our sum rule prediction fl3"5| ) can be tested experimentally, 
by a measurement of the form factor ratio R2 in B — ► D*£u transitions. Using 77(7/) = 0.6 
in flTT| ) we obtain 
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In Table | we compare this result to the predictions of some popular quark models, as well 
as to a recent QCD sum rule calculation of the weak decay form factors in the full theory. 
These models give values for R 2 which are substantially larger than ours. In particular, we 
note that at q 2 = 0, corresponding to the maximal velocity transfer, the models give Ri > 1, 
whereas we find R2 ~ 0.84. This discrepancy should not be too surprising. Since we have 
worked very hard to understand the origin of the symmetry-breaking corrections, we can 
hope that our refined sum rule analysis accounts for such effects in a much more detailed 
way than the naive quark models can. 

We end this paper with an interesting speculation. Although there is no reason to 
believe that it makes any sense to apply the heavy quark expansion to the D — > K*£ v decay 
amplitude, we might still believe in a "continuity of signs" and guess that the tendency 
R2 < 1 should persist, and most likely even become more pronounced, when we imagine 
changing the heavy quark masses from m;, and m c to m c and m s . This tendency is in 
fact very consistent with the experimental value of the form factor ratio obtained from an 
analysis of the joint angular distribution in D — > K*i v decays. Taking the weighted average 
of the results reported by the experiments E691 and E653 |35|], we get R% K * (q 2 = 0) = 
0.67 ± 0.23. Although we have no right to extrapolate fllT]) down to the strange quark 
mass, we take this observation as a confirmation of our prediction that symmetry-breaking 
corrections suppress i? 2 -Q 

In conclusion, we have presented the complete next-to-leading order QCD sum rule 
analysis of the subleading Isgur-Wise functions ^(v -v', /x) and T](v-v'), including in particular 
the two-loop perturbative corrections. We find that effects of order a s are very important 
and enhance the form factors. Over the kinematic region accessible in semileptonic decays, 
the renormalization-group invariant ratio rj(v ■ v') turns out to be essentially constant and 
equals 0.6 ± 0.2. This leads to an almost complete cancellation of the leading symmetry- 
breaking corrections to the B — > D iv decay rate at zero recoil, allowing for a reliable 
determination of | V c b\ from this decay mode. 
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x The authors of Ref. || apply HQET to D -> K^lv decays and obtain 77(1) 0.3 ± 0.4 (for 
A = 0.4 GeV) from an overall fit to the data. 
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TABLES 





our results 


ISGW §1] 


BSW H| 


KS|| 


sum rules [33| 


#2(1) 


0.80 


0.91 


0.85 


1.09 


0.95 


-^(l/max) 


0.84 


1.14 


1.06 


1.00 


1.05 



TABLE I. Bredictions for the form factor ratio i?2- The zero recoil limit y = 1 corresponds to 
q 2 = (tub — m£>*) 2 , whereas y max ~ L5 corresponds to q 2 = 0. 
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FIGURES 



FIG. 1. One- and two- loop perturbative contributions to the sum rule for the universal form 
factor £s(v ■ v',fj,). Heavy quark propagators are drawn as double lines. The wavy line represents 
the velocity-changing heavy quark current hi T iD a h. 

FIG. 2. Numerical evaluation of the sum rule (|30|): (a) dependence of the zero recoil form 
factor 77(1) on the Borel parameter for different values of the continuum threshold; (b) the function 
r)(y) for 0.6 < T < 1.2 GeV and 1.7 < lu < 2.3 GeV. 
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